OU-TAP 114 
UAB-FT 484 

The use of new coordinates for the template space in hierarchical search for 
gravitational waves from inspiraling binaries 

Takahiro Tanaka a,fe '* and Hideyuki Tagoshi™'' 

a Department of Earth and Space Science, Graduate School of Science 

Osaka University, Toyonaka 560-0043, Japan 
b IFAE, Departament de Fisica, Universitat Autdnoma de Barcelona 
08193 Bellaterra {Barcelona) , Spain 
*ta ma@vega.ess.sci.osaka-u.ac.jp 
^tagosh i@vega.ess.sci.osaka-u.ac.jp 

We discuss a method to analyze data from interferometric gravitational wave detectors focusing on the 
technique of hierarchical search to detect gravitational waves from inspiraling binaries. For this purpose, we 
propose new coordinates to parameterize the template space. Using our new coordinates, we develop several new 
techniques for two step search, which would reduce the computation cost by large amount. These techniques 
become more important when we need to implement a x 2 -test as a detection criterion. 



I. INTRODUCTION 

The interferometric gravitational wave detectors, such 
as, LIGO, VERGO, GEO600 and TAMA300 0, are now 
under construction. Especially, TAMA300 has already 
done the first large scale data acquisition for three days 
in September 1999 jy. Primary targets of these detectors 
are inspiraling binary neutron stars or black holes. These 
compact binaries can be produced as a consequence of 
normal stellar evolution in binaries. It is also suggested 
that they might also have been produced in the early 
universe. The analysis of the first 2.1yr of photometry 
of 8.5 million stars in the Large Magellanic Cloud by 
the MACHO collaboration suggests that 0.62to;2 of halo 
consists of MACHOs of mass O.5l[!jM0 in the standard 
spherical flat rotation halo model p|. If these MACHOs 
are black holes, it is reasonable to consider that they were 
produced in the early universe, and some of them are in 
binaries which coalesce due to the gravitational radiation 
reaction ||. Thus, it is expected that the observation of 
gravitational waves gives a definite answer to the question 
whether these MACHOs consist of primordial black holes 
or not. 

To search for gravitational waves emitted by these bi- 
naries, the technique of matched filtering is considered 
to be the best method. In this method, detection of sig- 
nals and estimation of binary's parameters are done by 
taking the cross-correlation between observed data and 
predicted wave forms. For this purpose, we need to pre- 
pare theoretically predicted wave forms, often called tem- 
plates. Generally, such templates depend on binary's pa- 
rameters such as mass, spin, coalescing time, phase, and 
so on. Since these parameters are continuous, what we 
really have to do is to prepare a template bank which 
consists of a finite number of representative templates. 

In this paper, we propose a new parameterization of 



two mass parameters of binaries. We show that the use of 
them has various advantages in performing the matched 
filtering. These parameters define two dimensional co- 
ordinates on the parameter space of templates. We can 
introduce a distance between two templates by using the 
cross-correlation between them. This distance defines a 
metric on the template space Q . We shall show that the 
metric in terms of our new parameters approximately 
becomes a flat Euclidean metric. Thus, it becomes very 
simple in these coordinates to determine the grid points 
corresponding to the bank of templates. The method 
how to construct the new coordinates is explained in the 
succeeding section. 

Requiring that the grid in the template space is suf- 
ficiently fine so as not to lose real events, the number 
of templates to be searched tends to be very large. Es- 
pecially, if we lower the minimum mass of binaries to 
be searched m m i n , the number of templates increases as 
m mili @ ■ When we search for gravitational waves from 
binary black hole MACHOs, we need to choose m m ; m 
sufficiently lower than the predicted mass of MACHOs 
~ 0.5M Q . For example, let us consider that the grid 
is chosen so that the correlation between nearest neigh- 
boring templates becomes 0.97. Then, in order to search 
for binaries composed of compact stars in the mass range 
between 0.2M Q and 1OM , the necessary number of tem- 
plates becomes ~ 2 x 10 5 for the "TAMA noise curve 
|| . The matched filtering with the sampling rate of 
3000Hz requires the data processing speed faster than 
80G FLOPS (FLoating Operations Per Second) for the 
on-line analysis. Now, such a powerful computation en- 
vironment may be available. However it is still very ex- 
pensive. Furthermore, there are various factors in real 
data analysis which increase the computation cost. One 
of them is the non-Gaussian nature of the detector noise, 
which we shall discuss in this paper. Hence, the com- 



putation cost can be much larger than that estimated in 
an ideal situation. Thus, it is required to develop some 
methods to reduce the computation cost. 

The technique of hierarchical search is thought to be a 
promising way to realize such reduction in the computa- 
tion cost ||. However, when we apply this technique to 
real data, the non-Gaussian nature of the detector noise 
mentioned above causes a trouble. As we shall explain 
later, a simple hierarchical search scheme does not work 
in the presence of non-Gaussian noise. To solve this diffi- 
culty, we propose some new computation techniques sup- 
plementary to the technique of hierarchical search. These 
new techniques depend very much on our choice of new 
mass parameters. 

This paper is organized as follows. In Section 2 we ex- 
plain the definition of our new coordinates which param- 
eterize the post-Newtonian templates. We also explain 
that the computation cost in the template generation 
process can be reduced by using our new coordinates. 
In Section 3 we discuss a difficulty in the hierarchical 
search, which has not been pointed out in literature, and 
explain a method to overcome this difficulty. Section 4 is 
devoted to summary and discussion. 

Throughout this paper, we use units such that New- 
ton's gravitational constant and the speed of light are 
equal to unity. The Fourier transform of a function h(t) 
is denoted by h(f), i.e., 
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We adopt the templates calculated by using the post- 
Newtonian approximation of general relativity @|. We 
use a simplified version of the post-Newtonian templates 
in which the phase evolution is calculated to 2.5 post- 
Newtonian order, but the amplitude evolution contains 
only the lowest Newtonian quadrupole contribution. We 
also use the stationary phase approximation,® whose va- 
lidity has already been confirmed in Ref. || . 

We denote the parameters distinguishing different tem- 
plates by M M . They consist of the coalescence time 
t c , the total mass m tot (:= m i + m 2)> the mass ratio 
?/(:= TOima/m^j), and spin parameters. The templates 
corresponding to a given set of M M are represented in 
Fourier space by two independent templates h c and h s 
as 



h = h c cos 4>q + h s sin cf>o , 
where 4>$ is the phase of wave, and 
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for < / < f max {M»), 
h c (M»J) = h a (M> t ,f) = 0, 

for f>fmax(M»). (2.4) 

Here J\f is a normalization constant, and 

Mf) = T t 6 i (MnCi(f), (2.5) 
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with 



II. NEW COORDINATES FOR TEMPLATE 
SPACE 

A. The noise spectrum and templates 

We assume that the time-sequential data of the detec- 
tor output s(t) consists of a signal plus noise n(t). We 
also assume that the wave form of the signal is predicted 
theoretically with sufficiently good accuracy. Hence the 
signal is supposed to be identical to one of templates ex- 
cept for the normalization of its amplitude. 

To characterize the detector noise, we define one-sided 
power spectrum density S n (f) by 



S n (f) = 2 (n(t)n(t + T))e 2mfT dT, 



(2.1) 



where ( ) represents the operation of taking the sta- 
tistical average. For the purpose of the present paper, 
the overall amplitude of S n is irrelevant. We adopt the 
"TAMA phase II" noise spectrum as a model, which is 
given by Q 
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Ci(/) = r 5/3 , c 2 (/) = r 1 ! c 3 (/) = r 2/3 , 

C 4 (/) = r 1/3 , Cs(/)=ln/. (2.6) 

We have quoted the expression for the case in which the 
spin vector of each star is aligned or anti-aligned with 
the axis of the orbital angular momentum. The spin 
parameters \s and \ a are related to the angular momenta 
of constituent stars S\ and Si by 
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where a plus (minus) sign is assigned to the angular mo- 
mentum when the spin is aligned (anti-aligned) with the 
orbital angular momentum. In (2.6), we have neglected 
the spin effects at 2.5PN order. Negative frequency com- 
ponents are given by the reality condition of h(t) as 



H-f) = UfY 



where * means the operation of taking the complex con- 
jugate. 

When we consider rather massive binaries, fmax must 
be chosen at the frequency below which the post- 
Newtonian templates are valid. On the other hand, 
when we consider less massive binaries, the maximum 
frequency f ma x is determined by the noise curve alone. 
In this case, we need to choose f max such that the loss 
of the signal-to-noise due to the discreteness of the time 
step, At c = 1/(2 fmax), is negligibly small. 



B. Template space in matched Altering and new 
parameters 

Here, we define the inner product between two real 
functions a{t) and b{t) as || 
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In the matched filtering, we define the filtered signal- 
to-noise ratio after maximization over 0o as 



p= \/(s,h c ) 2 + (s,h s y 



(2.11) 



We choose the normalization constant N to satisfy 
{h c , h c ) = {h a , h s ) = 1. 

Since the best fit value for parameters are not known in 
advance, we must filter the data through many templates 
at different points in the parameter space. In order to 
determine representative points in the parameter space, 
we have to know how much the value of p is reduced by 
using a template with different mass parameters from the 
best ones. Here, we adopt geometrical description of the 
template space Qj , and investigate which coordinates we 
should choose to simplify the strategy for determining 
representative points in the parameter space. 

In the following, we assume that the maximum fre- 
quency fmax is determined by the noise curve alone in- 
dependently of the template parameters. The effect of 
parameter dependence of f m ax is discussed at the end 
of this section. Although 9 % introduced above are func- 
tions of M M , we regard them as independent variables for 
a while. Hence, we parameterize templates like h c {9'), 
where we defined 0' as the set of parameters {t c , 9) set- 
ting 0' Q — t c and 9\ = 6i for i = 1, • • • , 5. 

The correlation between two nearby templates with 
different 9' is evaluated as 

{h c {9' + A6'),h c {9')) 
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where Aipg(f) — ipe+Ae{f) — ipe{f)- In the same manner, 
we have 
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Here it should be emphasized that these correlations de- 
pend on 9' only through AO' . 

Let us define new functions Q'(A9') and Q{A9) by 



g'{A9') := \{h c {9' + A9')M9')f 

+{h c {9' + A9'),h s {d')) 2 ] 1/2 , 



(2.14) 



and Q{A6) := maxAt c Q'{A9'), respectively. This func- 
tion Q is known as the match. We expand Q'{A9') with 
respect to A9' as 
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where we introduced a notation 
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We define a matrix G' by 

G'(AO') = 1 - G'^A&^AO'" 
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By definition, this matrix G' is a constant matrix inde- 
pendent of 9' , and it is determined once the noise spec- 
trum is specified. In order to take maximum of Q' with 
respect At c , we project G' on to the space orthogonal 
to t r as 



(2.18) 



The matrix G^- can be considered as a five dimensional 
metric analogous to the two dimensional one introduced 
in Ref. f§. 

Next, we orthonormalize as 



C — \ \ P a P a 



(2.19) 



where A Q are the eigenvalues and P a i is an orthogonal 
matrix composed of the eigen vectors of Gy . Rotating 
the axis further by using another orthogonal matrix Q A a , 
we define new coordinates of the five dimensional tem- 
plate space by 
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EE^ A « /2p "^> (^ = V-'.5). (2.20) 



Let us denote our five dimensional template space as 
r. In this paper, we assume that we can neglect the effect 
of spins of binary stars. Therefore, the actual template 
space to be searched becomes a two dimensional hyper- 
surface S in r. Sin ce 9 l are functions of m to t and m we 
find that Eq. fl2.20p defines a map from the actual tem- 
plate space parameterized by (mtot,v) to T. Then, this 
map naturally specifies S. 

One of the most important points that we wish to em- 
phasize in this paper is that the geometry of this two 
dimensional surface S becomes almost flat. Because of 
this fact, we can choose Q A a so that the x 1 and x 2 axes 
lie approximately on S. Taking into account the exten- 
sion of the area to be searched on S, we choose Q A a by 
solving the following set of equations, 



= (an, 0,0, 0,0), 

= (a 2 i,a 22 , 0,0,0), (2.21) 

where m m i n and m max are the minimum and the max- 
imum mass of templates, respectively. The directions of 
the other axes are not important here. (Hence, we do not 
specify how to choose them explicitly. ) 

By solving X 1 = x 1 (m to t,r]) and X 2 = x 2 {ra toU r\) 
for mtot and r], we obtain inverse functions mt t(X , X 2 ) 
and ^(X 1 , X 2 ), and we can use X 1 and X 2 as parame- 
ters for the template space instead of mtot and 77. Fur- 
thermore, we can define a map from (X 1 ^ 2 ) to T by 
x A = x A (m tot (X\X 2 ) : T 1 {X 1 .X 2 )) for A=3,4,5. 

In the following, we verify that the two dimensional 
hypersurface S is almost flat. First we check that the 



parameters x 



and x are approximately zero on any 



points on S. To show this, we plot the value of Sx :— 
\/ (x 3 ) 2 + (x 4 ) 2 + (a; 5 ) 2 as a function of m\ and m 2 in 
Fig.l. Here we set m,„j„ = 0.2Af Q and m max — 3M Q . 
We find that Sx are very small. This indicates that the 
surface S is almost flat and (X 1 ^ 2 ) can be regarded as 
Cartesian coordinates on S in a good approximation. 

The metric associated with the new coordinates is de- 
fined by 



9U :=E G 



dff 1 d6 j 
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Sab Y, \ dM » dx i 1 

A,B=1 |U=1 v 7 i>=1 

(J, J =1,2). 
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(2.22) 




FIG. 1. The value of Sx as a function of mi and m,2 in 
solar mass unit. 




To show the constancy of this metric, we plot the resid- 
uals gu - 81 j in Fig. 2 (a)- (c). 



FIG. 2. The values of (a)jn — 1, (h)g-22 — 1, (c)<?i2, as 
functions of mi and m,2 in solar mass unit. 
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These facts suggest the usefulness of the new coordi- 
nates. First of all, the flatness of the metric allows us to 
use a uniform square grid to generate the template bank. 
Besides, there are several advantages. As long as we con- 
sider a small area in T, x 3 ,x 4 and x 5 can be treated 
as constants. We can make use of this fact to develop 
an efficient algorithm to generate templates in frequency 
domain as we shall see in the succeeding subsection. 



C. An efficient algorithm to generate templates 

We can express 6 l as linear functions of x A by solving 
Eq.( [2.20 ) inversely. Thus, the phase function tp$(f) is 



also a linear function of x . Since the effect of variation 
off 

30 ^ 30 ^ 30 ) within a small area is negligibly small, the 
difference of the phase function 

AV>(/) = ^(xi + ajci ; xhax2)(/) ~ ip(x\x^){f ) (2.23) 

is almost independent of (X^X 2 ). Therefore, we can 
prepare the phase difference e lA ^^ for various values of 
(AX 1 , AX 2 ) in advance. Then we can calculate the tem- 
plate ft(xi+Axi,P + AX2)(/) just by multiplying the cor- 
responding phase difference by the template at (X , X 2 ) 
as 

^+axi,M»)(/) - e* AW) - (2.24) 

Hence, once we calculate one template, we do not have 
to call subroutines of the sinusoidal functions to gener- 
ate neighboring templates of it. Since the computation of 
sinusoidal functions is slow in many compiler, this algo- 
rithm significantly reduces the computation cost to gen- 
erate templates. 

Before closing this secton, we remark on the choice of 
the maximum frequency fmax- So far, we have been ne- 
glecting the fact that, in general, the maximum frequency 
fmax depends on the template parameters. When we 
consider binaries with large mass, the frequency at the 
last stable circular orbit becomes lower than the maxi- 
mum frequency which is determined by the shape of the 
noise power spectrum. Since our post-Newtonian tem- 
plates are no longer valid beyond the last stable orbit, 
the maximum frequency should be chosen below the cor- 
responding frequency. Even in that case, we think that 
our new coordinates are still useful by the following rea- 
son. The match determined with larger f ma x is likely to 
underestimate the correct value. Hence, the grid spac- 
ing determined by using our new coordinates tends to 
be closer than that determined by more accurate estima- 
tion. Therefore, to adopt constant f m ax in determining 
the template bank would be safe in the sense that it is 
less likely to miss detectable events. Although the num- 
ber of templates increases with the choice of constant 
fmax, such effect is very small. This is because the num- 
ber of templates with relatively large mass is not very 
large. Recall that the number of template are dominated 
by templates with small mass m ~ m m i n @ ■ 



III. NEW FAST ALGORITHM FOR 
HIERARCHICAL SEARCH 

If we try to search gravitational waves from coalescing 
binaries with mass M > 0.2Mq, we have to calculate 
the correlations (s(f), h(f)) for severalx 10000 templates 
with different mass parameters [To| ] R The computation 
cost to evaluate such a large number of correlations is 
very expensive. One promising idea to reduce the com- 
putation cost is the technique of hierarchical search. The 
basic idea of hierarchical search is as follows. At the first 
step, we examine the correlations with a smaller number 
of templates located more sparsely. In order not to lose 
the candidates of events, we set a sufficiently low thresh- 
old of the filtered signal-to-noise p at the first step. If a 
set of parameters (X 1 , X 2 , t c ) is selected as a candidate, 
we examine the correlations between the data and the 
neighboring templates by using a finer mesh. 

A simple scheme for two step search has been already 
discussed in Ref. ||. However, there seems to be a prob- 
lem in realizing the basic idea mentioned above. It has 
been pointed out that the distribution of the amplitude 
of the detector noise will not follow the simple stationary 
Gaussian statistics pf| . The non-stationary and non- 
Gaussian nature of the detector noise will produce a lot 
of events with large value of p. Namely, there seems to 
exist a non-Gaussian tail part in the distribution of p. In 
order to identify real events, we need to keep the expected 
number of fake events small by choosing the threshold of 
p as being sufficiently large. Hence, the existence of the 
non-Gaussian tail requires a larger value of threshold of 
p. This leads to a significant loss of detector sensitivity. 
To avoid such a loss, it was proposed to use a x 2 -test as 
a supplementary criterion. 

Here, x 2 is defined as follows. First we divide each 
template into mutual independent n pieces, 



and we calculate 

= ( g > ~ h fls)l = 4) x (*< k*>))> ( 3 - 2 ) 

with 

Then x 2 is defined by 
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1 If we take into account the effect of spins of binary stars, 
this number will increase about 3 times or more. We will 
discuss this issue in a separate paper Q . 
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This quantity must satisfy the x 2 -statistics with 2n — 2 
degrees of freedom and must be independent of p = 
long as the detector noise is Gaussian. 

However, as reported in Ref. pj| , events with large x 2 
in reality occur more often than in the case of Gaussian 
noise. There is a strong tendency that events with large 
X 2 have a large value of p on average. Thus, by changing 
the threshold of p depending on the value of \ 2 % w e can 
reduce the number of fake events without any significant 
loss of detector sensitivity. Hence, it will be necessary to 
implement the x 2 -test even in the simple one step search 
case. However, if we try to evaluate x 2 , the computation 
cost becomes more expensive^]. Thus, it is strongly de- 
sired to implement an efficient algorithm to reduce the 
computation cost. 

Now, we discuss a method of two step hierarchical 
search taking into account the presence of non-Gaussian 
noise. At the first step search, a large number of candi- 
dates for the second step search with large p value will 
appear due to the non-Gaussianity of noise. As is men- 
tioned above, in order not to lose the detector sensitiv- 
ity, it is desired to introduce the x 2 -test, and to keep the 
threshold of p small. Furthermore, by introducing the 
X 2 -test at the first step, we can reduce the number of can- 
didates for the second step search. Thus, the x 2 -test is 
also effective to reduce the computation cost for the sec- 
ond step. However, the x 2 -test at the first step increases 
the computation cost for the first step. This increase in 
the computation cost can be very large in the presence 
of non-Gaussian noise because the number of fake events 
which exceed the threshold of p at the first step is much 
larger than that expected in the case of Gaussian noise. 
Then, we must compute a lot of x 2 values at the first 
step. When we take into account these effects, the ad- 
vantage of the two step search, which is estimated to be 
about factor 30 in comparison with the simple one step 
search in the case of Gaussian noise 0, will be signif- 
icantly reduced or will be totally lost. Hence, in order 
to make use of the potential advantage of the two step 
search, we need other ideas to reduce the computation 
cost further. Here we present two new ideas of this kind. 



2 If we try to evaluate \ 2 naively, the computation cost nec- 
essary for the second step search simply becomes n times 
larger. Since n will be chosen as being O(10), the increase 
of the computation cost is unacceptable. If one calculates the 
values of \ 2 onr y f° r a few varieties of coalescence time at 
which a large value of p is achieved, the computation cost for 
X 2 might be kept small. In this case, we can use the direct 
summation instead of FFT to calculate the values of \ 2 ■ But, 
the question is for how many varieties of coalescence time we 
must calculate the values of \ 2 n °t to lose real events. This 
is not a simple question. If this number is sufficiently small, 
this naive strategy will work in the case of one step search. 
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TABLE 1 . Tables of maximum correlations for various choices 
of AX (a) with 5000Hz sampling and (b) with 1250Hz sam- 
pling. 

The first one is very simple. At the first step search, we 
can reduce the sampling rate of the data. The low sam- 
pling rate results in the reduction in the filtered signal-to- 
noise mainly due to the mismatch in t c . However such re- 
duction can be compensated by a very small change of the 
threshold of p. In the case of the "TAMA noise curve" , 
we can allow the sampling rate as low as about 1000Hz. 
The values of match between two templates with various 
AX are shown in Table. 1(a), where we adopted 5000Hz 
as the sampling rate. The same quantities for 1250Hz 
sampling are shown in Table. 1(b). Here, one of the tem- 
plates is considered as a normalized signal without noise, 
and the other as a search template. The signal is normal- 
ized to satisfy (h, h) = 1 with f max = 2500-ff z in both 
cases of Table. 1(a) and Table. 1(b). On the other hand, 
the search template is normalized with f max = 2500-ff z 
in the case of Table. 1(a) and with f ma x — 625Hz in the 
case of Table. 1(b). We find that the difference of the 
values of match between these two cases are very small 
especially for large |AX|. Therefore, detection probabil- 
ity for a fixed threshold of p is not significantly lost even 
if we adopt a rather small sampling rate at the first step 
search, at which we adopt a relatively large spacing for 
the template bank. The reduction in the sampling rate 
directly reduces the computation cost. The usual FFT 
routine requires effective floating point operations pro- 
portional to ~ N In N to compute the Fourier transform 
of the data with length N. Furthermore, for most of 
FFT routines and computer environments, the effective 
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FLOPS value for FFT is larger for smaller N. Thus, the 
reduction factor for the computation cost due to adopt- 
ing a smaller FFT length is much larger than one expects 
naively. 

The second idea is more important. What we need to 
evaluate is the correlation 



Z = (h, s) 



df 



fr, 



M/)S*(/) 



SJf) 
N r a*(/ J -)c i * x W*) 

fl /6 S n (fj) 



2irifjt a 



(3.5) 



where we take it into account that in reality we deal 
with a discrete time sequence of data with length N. 
A/ is given by the sampling rate divided by N, and 
fj := Af(j — N/2). The correlations for various values 
of t c are calculated simultaneously by taking the Fourier 
transform of the array defined by the quantity inside the 
square brackets in the last line of the above equation. 
This scheme is efficient enough when we do not have any 
guess about t c . However, when we perform the second 
step search we have a good estimate of t c at which the 
maximum correlation is expected to be achieved. We de- 
note it by t c . In this case, we need to evaluate Z only 
for t c which are close to t c . Also for templates, we have 
a good guess for the mass parameters, X. Thus, we need 
to calculate the correlation Z only for a cluster of the 
templates neighboring to X. 

Once t c and X are specified, we can rewrite the above 
expression in a very suggestive form as 



N/m-l m(fe+l/2) 

z^MAf e E A i xB i> 

k=l j=m(k-l/2) + l 



(3.6) 



with 



B . _ e iAV(/ J )+2 7 ri/ J At c 



(3.7) 



tr 



tr 



Here Aip(fj) := ^ x (/,) - Vx(/i) and At < 
We have also introduced to as a certain integer which 
divides N. As long as both |X — X| and \t c — t c \ are 
sufficiently small, the factor B is a slowly changing func- 
tion of frequency. Hence, unless to is not larg e, B can be 
moved outside the second summation in Eq.(3.6). Then, 
introducing 



m(ft+l/2) 

= E 

j=m(fc— l/2)+l 



(3.8) 



we obtain 



N/m-l 

Z » ALA/ A 'k x B mk 

k=l 



N/m-l 

-NAf ]T 



fe=i 



[A' k e' 



x e 



2irif' k At c 



(3.9) 



where 



iAf(k 



N 

2m I 



The expression in the last 



line can be evaluated by applying the FFT routine to the 
array defined by the quantity inside the square brackets. 
The correlation between two templates for various values 
of AX and At c are calculated by using this method. The 
results are shown in Table. 2 for N/m = 1024, 2048 and 
4096. We find that N/m can be taken as small as 2048 
without significant loss in accuracy. 



(a) 



Ai e (sec) 


m = 4096 


m = 2048 


m = 1028 


0.0000 


1.000 


1.000 


1.000 


0.0128 


1.000 


0.999 


0.994 


0.0248 


0.999 


0.994 


0.975 



(b) 



Ai e (sec) 


m = 4096 


m = 2048 


m = 1028 


0.0000 


0.765 


0.765 


0.765 


0.0128 


0.764 


0.763 


0.760 


0.0248 


0.763 


0.760 


0.746 



(c) 



Ai e (sec) 


to = 4096 


m = 2048 


m = 1028 


0.0000 


0.774 


0.774 


0.774 


0.0128 


0.774 


0.773 


0.769 


0.0248 


0.773 


0.769 


0.755 



TABLE 2. Tables of maximum correlations for various choices 
of m and At e = t e -t ( e c) with (a) AX = (0,0), (b)AX = (1,0) 
and (c)AX = (0,2). 

Furthermore, as an advantage of our new coordinates, 
the factor e l ^^ mh ' > can be well approximated by the 
one obtained by setting Ax 3 — Ax 4 — Ax 5 = 0. It 
means that this factor is almost independent of the val- 
ues of (A 1 , A 2 ). Thus, we have to calculate this factor 
only once at the beginning of the second step search. 
Since the array A' k is independen t of AX, the quantity 
inside the square brackets in Eq. (3.9) for various values 
of AX is simply given by A' k times the pre-calculated fac- 
tor e lA ^^ mk \ This fact manifestly leads to an additional 
reduction in the computation cost. 

The same technique can be used to evaluate \ 2 : i- e -; 
Z^\ just by replacing the array A' k with the same quan- 
tities multiplied by an appropriate window function. 



IV. CONCLUSION 

We discussed a method of analyzing data from inter- 
ferometric gravitational wave detectors to detect gravita- 
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tional waves from inspiraling compact binaries based on 
the technique of matched filtering. We described a brief 
sketch of several new techniques which would be useful 
in hierarchical search of gravitational waves. 

First, we proposed new parameters which label tem- 
plates of gravitational waves from inspiraling binaries. 
These new parameters are chosen so that the metric on 
the template space becomes almost constant. We found 
that the template space can be well approximated by two 
dimensional flat Euclidean metric. Thus, by using these 
parameters as coordinates for the template space, the 
problem of the template placement becomes very sim- 
ple. Say, we can use a simple uniform square grid to 
specify the grid points for the bank of templates. Fur- 
thermore, we found that, by using new parameters, we 
can introduce an efficient method to generate templates 
in frequency domain. The reduction in the computa- 
tion cost is achieved by using the property of our new 
coordinates that one template can be translated into an- 
other with different mass parameters by just multiply- 
ing pre-calculated coefficients to the original template. 
Therefore, we can generate a set of templates from one 
template avoiding calculation of the sinusoidal functions. 

Next, we discussed a method of two step hierarchical 
search. Due to the non-stationary and non-Gaussian na- 
ture of the detector noise, we will have to introduce a 
X 2 -test when we analyze real data. When we take into 
account this fact, it becomes very difficult to obtain large 
reduction in the computation cost by applying naive two- 
step hierarchical search strategy. To solve this difficulty, 
we proposed two new techniques to reduce the computa- 
tion cost in the two step search. One is to use a lower 
sampling rate for the first step search. By using this tech- 
nique, we can reduce the length of FFT by factor 2 or 4 
keeping the loss of correlation within an acceptable level. 
The second technique, which is more important, makes 
use of the fact that a good guess for the coalescence time 
and the mass parameters has been obtained as a result 
of the first step search at the time when we perform the 
second step search. We found that the length of FFT 
for the second step search can be reduced down to about 
2048. 

Based on the new techniques discussed in this paper, 
we have developed a hierarchical search code to analyze 
data from the TAMA300 detector. The details of this 
code and the result of the analysis of the first TAMA300 
data will be presented elsewhere. 
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Some of the numerical calculation were done by using the 
gravitational wave data analysis library GRASP [pL3|| . 
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